In [7]:
%pylab inline
import phreeqpython
pp = phreeqpython.PhreeqPython('phreeqc.dat')
In [23]:
print("Lower limit", pp.add_solution({'Al': '0.2 mg/kgw', 'pH': '4 Gibbsite'}).pH)
print("Upper limit", pp.add_solution({'Al': '0.2 mg/kgw', 'pH': '8 Gibbsite'}).pH)